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Over the past few years bulk pyrochlore iridates of the form A 2 lr 207 (where A is a rare earth 
element, Ir is iridium, and O is oxygen) have been studied as model systems for investigating the 
interplay of electronic correlations and strong spin-orbit coupling, particularly with the aim of finding 
correlation-driven topological phases. In this work, we use cellular dynamical mean field theory 
(CDMFT) to study effects of electronic correlations beyond Hartree-Fock theory in thin films of 
pyrochlore irradiates grown along the [111] direction. We focus on the bilayer and trilayer systems, 
and compute the phase diagrams of these systems as a function of electron-electron interaction 
strength, which is modeled by an on-site Hubbard interaction. By evaluating the Z 2 invariant 
and Chern number using formulas based on the single-particle Green’s function and the quasi¬ 
particle effective Hamiltonian, we show that on-site correlations can drive an interaction-induced 
topological phase transition, turning a time-reversal invariant topological insulator and a nearly 
flat band metal to a correlated Chern insulator (Cl) in bilayer and trilayer systems, respectively. 

By comparing with the Hartree-Fock results, the CDMFT results show that quantum fluctuations 
enhance the robustness of the interaction-driven Cl phase in the thin films. Furthermore, our 
numerical analysis of the quasiparticle spectrum reveals that the topological phases we find in our 
many-body calculations are adiabatically connected to those in the single-particle picture. 

PACS numbers: 


I. INTRODUCTION 

Over the past decade, time-reversal invariant topo¬ 
logical insulators (TI) in two and three dimensions 
have received significant theoretical and experimental 
attention^-— Recently, significant effort has been made 
to investigate the role of electron-electron interactions in 
topological states of of matterThough it is gen¬ 
erally understood that topological phases described by 
electronic band structure are robust to weak electron- 
electron interactions, the full many-body problem of a 
strongly interacting system remains far from completely 
understood One possible effect of electronic interac¬ 
tions is magnetic order that carries with it a topological 
phase transition (either from a non-topological system 
to a topological one, or from a TI to a non-topological 
magnetically ordered state) 

Among the real-material proposals for systems ex¬ 
pected to exhibit topological properties, transition metal 
oxides (TMO) are a promising candidate for realizing 
nontrivial interacting topological phasesM^— Among the 
bulk (as apposed to thin-film) TMO that have been theo¬ 
retically studied in the context of topological phases, the 
irradiates have been particularly singled In 

these materials, the iridium 5 d electrons typically domi¬ 
nate the states near the Fermi energy. The d-orbitals are 
usually correlated and the large atomic number of irid¬ 
ium means that spin-orbit coupling can be significant as 
well. Thus, the iridates are a natural platform for inves¬ 
tigating the interplay of strong spin-orbit coupling and 
correlation effect s^^ 

In addition to the search for three-dimensional topo¬ 
logical phases in bulk TMO oxides, thin films (as thin 
as a bilayer) of TMO have also received significant at¬ 
tention recently;^ A number of these studies have fo¬ 


cused on the perovskite structure ABO3 (where A is 
a rare earth element, B is a transition metal, and O 
is oxygen) where thin films are grown along the [111] 
direction.— In this work, we are interested 
in thin films of the pyrochlore irididate ^2^207 (where 
A is a rare earth element, Ir is iridium, and O is oxy¬ 
gen) grown along the [ 111 ] direction 4 ^“— The thin film 
geometry of ^2^207 allows one to investigate the effects 
of reduced dimensionality on the interplay of correlations 
and strong spin-orbit coupling. Along the [111] direction 
(or equivalent directions), Ir^+ ions form alternating lay¬ 
ers of kagome and triangular lattices. See FigH] for a 
visualization of the geometry of the Ir atoms in the bi¬ 
layer and trilayer geometry. 

Our main objective in this work is to go beyond the 
Hartree-Fock (HF) approximation^!^ in the study of the 
effects of reduced dimensionality on the interplay of cor¬ 
relations and strong spin-orbit coupling in [111] grown 
pyrochlore iridate thin films. The intrinsic spin-orbit cou¬ 
pling of the Ir local moments generally breaks the spin- 
rotational invariance^i^ and therefore opens the possi¬ 
bility of long-range magnetic order at non-zero tempera¬ 
tures. Therefore, we do not expect qualitative changes to 
the magnetic order upon the inclusion of spatial fluctua¬ 
tions on top of the magnetic order predicted in HF the¬ 
ory. On the other hand, we wish to better understand 
how the quantum fluctuations in time may impact the 
HF predictions of interacting topological phases (namely 
the Chern insulator) in these systems. In particular, we 
would like to know whether the temporal fluctuations can 
be so severe that interaction-generated topological phases 
predicted by HF theory disappear. For Chern insulators 
induced “purely” by electron-electron interactions (that 
is, the underlying Hamiltonian has no intrinsic spin-orbit 
coupling at the non-interacting level), it appears that 
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temporal fluctuations can indeed destroy the topological 
phase^^ In our model-that does include a finite spin- 
orbit coupling at the non-interacting level-we reach the 
opposite conclusion: Temporal fluctuations do not gen¬ 
erally destroy an interaction-driven topological phase in 
two-dimensions. 

To investigate the physics of temporal fluctuations in 
our thin-film systems, we employ cellular dynamic mean- 
field theory (CDMFT) j^^i^^ We focus on two thin film 
systems of A2IT2O7 grown along the [111] direction: a 
bilayer consisting of a kagome and a triangular lattice 
(KT), and a trilayer consisting of triangular-kagome- 
triangular (TKT) lattices, as shown in Fig.[T][a) and (b). 



FIG. 1: (Color online) The lattice structure of the (a) bilayer 
(TK) and (b) trilayer (TKT) pyrochlore iridate thin films 
grown along the [111] direction. Each lattice site denotes an 
Ir^^ atom. The light green dots form a coplanar kagome lat¬ 
tice, whereas the pink dots form a coplanar triangular lattice, 
(c) shows the hexagonal Brillouin zone of the TK and TKT 
systems. 

Because our interest in this work is primarily the in¬ 
vestigation of temporal fluctuations on interaction-driven 
topological phases in two dimensional models (as opposed 
to a first-principles prediction of interaction effects in 
real material systems^) , we will simplify the problem by 
restricting the Hilbert space of our Hamiltonian to the 
= 1/2” model, appropriate for the limit of a large 
spin-orbit coupling^^ While the large spin-orbit coupling 
limit may not be reached in real materialswe expect 
the central physical results of our study to be unchanged 
upon the inclusion of the full t2g sub-space of the irid¬ 
ium bd orbitals. The reduced Hilbert space allows us to 
carry out the CDMFT study; the inclusion of the full t2g 
sub-space is beyond the numerical capabilities presently 
available. Because the bandwidth of the states near the 
Fermi energy is rather narrow, even a moderate interac¬ 
tion can put the system in a strongly correlated regime. 

Our CDMFT results show that the correlation effects 
can drive a topological phase transition from a TI to 
a magnetic Chern insulator (Cl) with nontrivial Chern 
number C = 1 in the bilayer (KT) system. In the tri¬ 
layer (TKT), we find that many-body interaction effects 
on the nearly flat band near the Fermi surface can in¬ 
duce a topological phase transition from a nonmagnetic 


conductor (C) to an interacting Cl with Chern number 
± 1 . Both thin film cases show that the ground state is a 
trivial magnetic insulator in the strong interaction limit. 
Compared to the previous HF results^ our results sug¬ 
gest that moderate quantum fluctuations (captured by 
CDMFT) can stabilize the interaction-driven Cl phases 
in both bilayer and trilayer systems, effectively enlarging 
the parameter space where they appear. 

Our paper is organized as follows. In SeclHl we first 
introduce a simplified model Hamiltonian for the thin 
films of pyrochlore iridates ^2^207. In Sec JHIl we briefly 
describe the CDMFT method, and explain how to eval¬ 
uate the Chern number and the Z2 invariant. In Sec. 
lYl we present our numerical results for the bilayer and 
trilayer cases with the exhaustive phase diagrams and 
corresponding magnetic configurations. Finally, we sum¬ 
marize our work in Sec. El 


II. MODEL HAMILTONIAN 

In the pyrochlore oxides H2H2O7, with A a rare earth 
element such as Y or La, the H-site has a vanishing 
magnetic moment and the physics is dominated by the 
transition metal B ions. In this work, we apply the 
model Hamiltonian proposed for the three-dimensional 
materials2^i^ to the quasi-two-dimensional thin film 
systemsEffectively integrating out the oxygen orbitals 
involved in indirect hopping processes between H-site 
ions, one obtains a model that only involves the H-sites. 
We will use such a model in this work. The bd atomic 
orbitals in Ir^+ are subject to a cubic crystal field, which 
splits the bd orbitals into Cg and t2g manifolds. The Cg 
manifold typically lies 2-3 eV above the t2g manifold^^ 
Spin-orbit coupling (SOC) further splits the t2g manifold 
into a Jeff = 1/2 doublet and a Jeff = 3/2 quadruplet. 
At infinite SOC strength, the Jeff = 1/2 states form a 
low energy manifold (the Jeff = 3/2 levels are inactive 
since they are fully occupied below the Fermi level) for 
the thin films. 

In the pyrochlore iridates ^2^207, four Ir^+ ions form 
a tetrahedron in the unit cell of a face-centered cubic Bra- 
vais lattice and the Jeff = 1/2 manifold is half-filledIn 
addition, we include an on-site Coulomb repulsion within 
the Jeff = 1/2 pseudospin space to obtain the effective 
Hamiltonian for our thin films 

H= E + (1) 

(Ri ,Rj) ,crcr' 

Ri,cr Ri 

where cr-ct annihilates an electron with pseudospin a at 
the ith site of the Bravais lattice vector R. The site in¬ 
dex i runs from 1 to 4 in the bilayer KT lattice [in Fig. 
dja)], while it runs from 1 to 5 in the trilayer TKT lat¬ 
tice within a single unit cell [in Fig. mb)]. The hopping 
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parameter Toxy arises from the oxygen-mediated hopping 
between nearest neighbor (NN) Ir^+ atoms with ampli¬ 
tude t oc where Vpd is the tunneling amplitude 

between p-orbitals of the oxygen and d orbitals from Ir; 
A is the energy difference between the two atomic orbital 
levels. The hopping parameter Tdir is the direct NN Ir-Ir 
hopping due to the direct overlap between the extended 
bd orbitals, which depends on to- and from a-ir bond¬ 
ing between the d-orbitals#^ The chemical potential ji is 
adjusted so that at half filling, each site is has one parti¬ 
cle on average in the Jeff = 1/2 bands. The values of the 
tight binding parameters in the matrices [Toxy] and [J^] 
have been explicitly expressed in Appendix A of Ref. [T 7 | . 
For simplicity, we choose hopping between oxygen p or¬ 
bitals and iridium d orbitals, tpd = 1, as the unit of 
energy, and set tyr = —2tcjlS for both bilayer and trilayer 
system in order to compare with bulk system o^^i^^ and 
thin film results 4 i 

Following earlier work on pyrochlore iridiates^i^ we 
select different values of the tight-binding hopping pa¬ 
rameter to- to explore different regimes of the parameter 
space potentially relevant to experimental systems. Be¬ 
cause of the underlying triangular lattice in the (111) 
planes, the Brillouin zones of both the KT and the TKT 
systems are hexagonal-see Fig. [He). 

For the bilayer system, we choose = 2 and present 
the noninteracting band structure and density of states 
(DOS) in Fig. [ 2 ][a). The KT lattice breaks inversion sym¬ 
metry and the noninteracting band structure has an en¬ 
ergy gap around the Fermi level at half filling. A Kramers 
degeneracy is found at time reversal invariant momenta 
(TRIM), i.e., F and M in Fig. [Il[c). By calculating the 
Z2 invariant using the Fu-Kane formula!^ or with Fukui’s 
method^ one finds the system is a topological band in¬ 
sulator for = 04 ^ For this value of ter, we can thus 
calculate within the CDMFT formalism how interactions 
drive the system away from a Z2 TI as the parameter U 
is increased. 

On the other hand, we choose = —1 for the trilayer 
case. In this parameter regime, the non-interacting band 
structure, depicted in Fig. Efb), shows nearly degenerate 
flat bands close to the Fermi level, such that at half filling 
it is metallic. As is revealed by previous HF results^ 
however, a small but finite interaction U will open a gap, 
and drive the system to a CL We wish to study the fate 
of this transition within the CDMFT formalism. 


III. CDMFT FORMALISM AND EVALUATION 
OF TOPOLOGICAL INVARIANTS 

CDMFT is an extension of the single-site dynamical 
mean-field theory method^ to include spatial correlation 
effects within a super-lattice unit cel L^^i^^ As a result, 
this method can capture the short-distance effect of quan¬ 
tum fluctuations neglected in single-particle approxima¬ 
tion and has already been successfully applied to study 
other complex oxide o^^i^^~ — and the three-dimensional 


(a)TI:t^=2.0 (b)C:t^=-1.0 




FIG. 2: (Color online) The noninteracting band structure 
along high symmetry directions in the Brillouin zone and den¬ 
sity of states (right side of each figure shown in red) for (a) 
bilayer (KT) at ta = 2; and (b) trilayer (TKT) at ta = — 1. 
At half filling, (a) is insulating whereas (b) is metallic. 

bulk (and slab) pyrochlore iridates^^S In general, the dy¬ 
namical mean-field theory method reduces lattice prob¬ 
lems with infinite degree of freedom to a type of An¬ 
derson impurity problem, in which a cluster of Nc sites 
hybridizes with electron bath sites (environment): 

Himp = ^ucf^vct' T ^ 

pva fj, 

+ E (W Acm. + 

plaa' 

T ^ ^ (2) 

l(T 

where the greek symbols //, z/ = 1, Ac label cluster sites; 
I = 1,..., A5 labels bath sites, and cr, a' are pseudo-spin 
labels. The electron operators (c^a) apply to the 

cluster sites, whereas aj^ {aio-) to the bath sites. The hop¬ 
ping integrals and chemical potential within the cluster 
are incorporated by the matrix elements of A, which are 
obtained from the tight-binding parameters [Toxy] and 
[Tdir]. The V’s and e’s are bath parameters describing 
hybridization between the clusters and bath sites, and 
on-site energy levels, respectively. The values of the V, 
the e and the lattice Green’s functions of Eq.([2]) are nu¬ 
merically determined via an exact diagonalization and 
self-consistency procedure. 

The self-consistency procedure is as follows: an initial 
input of the bath parameters V and e is given in Eq. (|2]) to 
solve the impurity Hamiltonian. Erom the cluster impu¬ 
rity Hamiltonian, we compute the cluster Green’s func¬ 
tion G{iuj) as a 2Ac x 2Ac (number of sites = Ac and 
number of pseudospin degrees of freedom = 2) matrix, as 
well as the cluster self-energy 

Ec(iw) = [g{ioj)]~^ - ( 3 ) 

where Q is the non-interacting cluster Green’s function. 
In CDMET, it serves as the Weiss field describing the 
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coupling of the cluster to the environment. Then the 
lattice Green’s function is coarse-grained as 

Gioc(iw) = ^ X] + M)i - ^(k) - Sc(«w) , (4) 

kSBZ 


where k is the momentum in the Brillouin zone, and t(k) 
is the Fourier transformed hopping integral. The self- 
consistent loop is closed by calculating the new Weiss 
function 


[Qnew{il^)] ^ = [Gloc(iw)] ^+Sc(iw), ( 5 ) 

and requiring [0new(^<^)] — fo within a prescribed 

accuracy. In updating the V and the e, we fit the new 
Weiss function on the Matsubara axis with conjugate gra¬ 
dient methods and chose an artificial low temperature 
ksT = 0.01 (in units of t) to simulate the zero tem¬ 
perature case. This value is much smaller than other 
characteristic energy scales in the problem, including the 
energy spacing due to the discretization of the momen¬ 
tum space in the Brillouin zone, and thus approximates 
the zero temperature limit well. 

In our work, we chose a tetrahedron as a unit cell clus¬ 
ter for the bilayer case, i.e. A^c = 4 ; while for the tri¬ 
layer case, the unit cell cluster is chosen as a tetrahe¬ 
dron with an external (pink) site linked to one of its 4 
corners, i.e. Nc = 5 [cf. Fig. [TJa) and (b)]. In order 
to solve Eq.([2]), exact diagonalization with the Lanczos 
algorithrr^ is employed to get the ground state proper¬ 
ties of the cluster, such as the magnetic order parameter 
and single-particle Green’s function. In the bilayer case, 
we compare = 4 and A^^ = 8 to determine how the 
phase diagram depends on the finite bath sites^; In the 
trilayer case, we choose A^^ = 5 . 

Next, we briefly describe how we numerically calculate 
topological invariants using the CDMFT method. In¬ 
terested readers can see more details in Appendix In 
a non-interacting topological band insulator, the Ghern 
number is equal to the TKNN number— 

^ = T f:. (S) 

Q;=filled 

where fij = dai (k) /dkj — daj (k) /dki is the Berry cur¬ 
vature, and ai(k) = —i{ua(k)\-^\ua(k)) is the Berry 
connection. The summation is over filled band indicies. 
Here, |ua(k)) denotes the a-th eigenstate of the nonin¬ 
teracting Hamiltonian H(k). There is no inversion sym¬ 
metry in the pyrochlore oxide [111] grown bilayer system. 
To calculate the Z2 index, (—1)^^^ generic way, we 
resort to the approach given by Fukui^ in the half-BZ, 

P2 = [ d?kf 12 — f dk • a mod 2. ( 7 ) 

Jbz JdBZ 

Note that in the Z2 formalism, a time-reversal smooth 
gauge has to be chosen as 


(8) 


where T is the time-reversal operator. 

For a general interacting system, one is no longer 
able to use the non-interacting invariants given in Eq.(|6]) 
and Eq.([ 7 j). Instead, one must express the invariant 
in terms of the single-particle Green’s function at zero 
frequency j^^i^^ With the CDMET method, we can easily 
obtain the zero-frequency Green’s function G{iuj = 0, k) 
at self-consistency. The zero-frequency Green’s func¬ 
tion allows one to define a “topological Hamiltonian”— 
from which one can compute the corresponding invariant, 

Ftop(k) = -G-i(0,k). (9) 

One diagonalizes the “topological Hamiltonian” to obtain 
the eigenvalues 


Htop\a,k) = ^ia{k)\a,k), ( 10 ) 

where a is a “band” index. Tilled “bands” are selected as 
eigenvalues with //^(k) < 0 . Then the Berry connections 
ai(k) are constructed by replacing the non-interacting 
Hamiltonian eigenstates |uQ,(k)) with the topological 
Hamiltonian eigenstates |<a,k), and plugging the result 
into Eq.® and Eq.© to evaluate the Ghern number 
and Z2 invariant, respectively. 

Eor weak to moderate strength interactions, quasipar¬ 
ticles exists with a finite lifetime. Because of the lifetime 
broadening of the states, it can be difficult to determine if 
a gap exists by directly evaluating the interacting quasi¬ 
particle spectral function A(cL;,k) = + i?7,k): 

The quasiparticle bands will have a finite width in en¬ 
ergy and are sensitive to the value of 77. In order to deal 
with this numerical issue, we determine the “effective” 
quasi-particle dispersion from the quasiparticle effective 
Hamiltonian Eq. (jB 14 p by expanding the local self-energy 
up to first order in frequency. The details are explained 
in Appendix [0 We emphasize that this effective Hamil¬ 
tonian mainly captures the quasi-particle features around 
the Eermi energy (cj = 0 ). 


IV. NUMERICAL RESULTS 

In this section, we present our GDMET results. In the 
HE theory, the single-particle band structure is renor¬ 
malized by the Hartree and Eock terms in the mean-field 
Hamiltonian^ and the self-energy of the single-particle 
Green’s function depends only on the momentum k (be¬ 
ing independent of the frequency uj). Therefore, we re¬ 
examine the fate of topological phases and transitions 
obtained by HE with CDMET. In the latter approach, 
quantum fluctuations and correlations are included non- 
perturbatively in the self-energy within the cluster. How¬ 
ever, in CDMET calculations, one usually obtains a hys¬ 
teresis behavior when the phase transition is of the first 
order. Namely, the critical value of U (say Ud) when the 
system undergoes a transition from phase A to phase B 
by increasing U does not coincide with the one (say Uc2) 
when the system undergoes a transition from phase B 


|u,(-k)) =T|u,(k)), 
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FIG. 3: (Color online) Bilayer (KT) phase diagram for to- = 2. 
(a) iV;, = 4 CDMFT results; (b) Nb = S CDMFT results; (c) 
HF results. In (a), the net magnetic moment defined in the 
unit cell m (black circles) and averaged magnetic moment per 
site ffi (blue triangles) are plotted as a function ofU. TI: time- 
reversal invariant topological insulator; MI: trivial magnetic 
insulator; Cl: topological Chern insulator. 


to phase A by decreasing U. In order to determine the 
complete phase diagram, we numerically calculated and 
compared the free energy per site when both magnetic 
and non-magnetic solutions coexist for a particular value 
of U. The phase is assigned to the state with the lower 
free energy. 


A. Bilayer system 


(a) Tl: t^=2.0,U=2.0 (b) Cl: t^=2.0,U=4.0 (c) MI: t^=2.0,U=6.0 
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FIG. 4: Spectral weights along high symmetry lines and local 
density of states for the bilayer (KT) with to — 2^ Nb — 8 . 
The broadening factor is 77 = 0.02. The yellow dots show 
quasiparticle spectrum obtained from Ea. (IB14D . (a) TI: U = 
2.0; (b) Cl: U = 3.0; (c) MI: U = 6.0 



remains in the time-reversal symmetric (TRS) TI phase. 
The spectral functions and the quasiparticle band struc¬ 
tures (see Appendix [B]) computed from = S CDMFT 
are plotted in Fig. lU^a), while A5 = 4 gives similar re¬ 
sults. We can see from Fig. SJ^a) that the Kramers de¬ 
generacy for quasiparticle bands remains intact and the 
finite band gap leaves the Z2 invariant z/ = 1 unchanged. 
In this case, there is no magnetization induced, and the 
phase is adiabatically connected to the noninteracting TI 
phase. Thus, the TI is stable against moderate electronic 
correlation, as would be expected for a system with a gap 
in the excitation spectrum about the Fermi energy. 

For U > 3.0 (in CDMFT), TRS is broken and the 
system becomes a magnetic insulator. At Uc ~ 3 . 0 , 
a magnetization jump appears, which indicates a first- 
order phase transition from a non-magnetic to a mag¬ 
netic insulator. The net magnetic moment, in units of 
gfiB depicted by black circles in Fig. [ 3 fa), is defined as 

i 


where i runs over all sites within a unit cell. On the other 
hand, in Fig. Ha), the averaged magnetic moment per 
site, defined as 


The A5 = 4 and = 8 CDMFT phase diagram to¬ 
gether with the magnetization in the unit cell are pre¬ 
sented as a function of U in Fig. [ 3 l^a) and (b), respec¬ 
tively. While our results are for the parameter to- = 2 , we 
also confirmed that a similar phase diagram is obtained 
with ter = 1 . For comparison, we also show the phase 
diagram obtained by Hartree-Fock theory Fig. ^c). 

Both the CDMFT and HF approaches show that there 
exist three different phases upon increasing U, and that 
the phase boundaries are quantitatively similar, except 
for the region of the Chern insulator (Cl) which is much 
larger when quantum fluctuations are included. There 
exists a finite range of values, 0 < U <3.5inA5 = 4 and 
0 < [/ < 3.0 in Nb = 8 (CDMFT), for which the system 


shows a monotonicaly increasing magnetization upon in¬ 
creasing U. A similar behavior is also captured using 
the HF theory shown in Fig. [ 3 fb). In this regime, the 
ground state is a trivial magnetic insulator (MI) due to 
the trivial Chern number (7 = 0 . 

Upon increasing (7, a finite range that harbors the 
interaction-induced Cl phases is observed at (7 ^ 3 . 8 — 4.8 
for = 4 , (7 - 3.8 - 4.7 for A5 = 8 by CDMFT and at 
U ~ 3.63 — 3.65 by HF theory. In comparison between 
A5 = 4 and = 8, it is found that the magnetic phase 
boundary is shifted by around 0.5 in U. We ascribe this 








































FIG. 5: Magnetic configurations of the bilayer system within 
a unit cell at = 2 and at (a) U = 4.0 (Cl) and (b) U = 5.0 
(MI). Both magnetic configurations resemble the antiferro¬ 
magnetic ordering discovered in Ref. [^1^ : where moment 
1 (2) intersects with moment 4 (3), or moment 1 (2) inter¬ 
sects with moment 3 (4). To compare with the bulk case,— 
(c) takes ta = —l,U = 4.0 and is a trivial MI. The configura¬ 
tion (c) is close to the “All-in/All-out” (AIAO) order obtained 
in the three-dimensional bulk,— but with a non-vanishing net 
moment pointing in the (111) direction. The AIAO order will 
have zero m but nonzero ffi. 


shift to the effect of finite bath sites. The spectral func¬ 
tions and the quasiparticle dispersion for the Cl from 
Ni) = S are plotted in Fig. IH^b). To avoid the broad¬ 
ening of the spectral function by the imaginary part of 
the self-energy, we identify the gap around Fermi level by 
examining the band gaps from both topological Hamilto¬ 
nian Eq.®, and the quasiparticle effective Hamiltonian 
EQ. (|B14p . In both themes, the band topologies are well- 
defined and a nontrivial Chern number C = +1 is found 
for half-filled bands. One can observe from Eig. |4] that 
the quasi-particle dispersion matches well with the ridges 
of spectral weights along high symmetry lines, according 
to our definition of quasi-particle effective Hamiltonian 
EQ. (jB14p . The Kramers degeneracy at the TRIM points 
is lifted in the presence of finite magnetic moments. More 
interestingly, the Cl phase “survives” for wider range of 
U values in the CDMET phase diagram than in the HE 
phase diagram. This observation suggests that quantum 
fluctuations stabilize the interaction-induced Cl phase. 

The magnetic configurations for the Cl and MI for 
Ni) = S are illustrated in Eiglbfa) and (b), respectively. 
Similar configurations are also found for = 4. Eor 
the parameters we have studied, the magnetic config¬ 
urations in the bilayer case are similar to that in the 
bulk pyrochlore oxidesdenoted as Fs^ except for 
the non-vanishing net magnetization in the tetrahedron 
(m 7 ^ 0 ) due to spatial anisotropy between the in-plane 
and out-of-plane directions. This pattern is called the 




FIG. 6: Trilayer (TKT) phase diagram for ter = —1. (a) 

GDMFT results; (b) HF results. In (a), the net magnetic 
moment defined in the unit cell m (black circles) and averaged 
magnetic moment fh (blue triangles) are plotted as a function 
of U. G: trivial nonmagnetic conductor; M: trivial magnetic 
conductor. 


antiferromagnetic ordering in the previous studies of the 
bulk system^^^i^ In the Cl phase of Fig. [3fb) (A^^ = 8 ), 
we found a discontinuity in the magnitude of net mag¬ 
netization around U ~ 4.0 due to a sudden drop of the 
ratio between the magnitude of moment 3 and moment 
1 in Fig. [5fa), while the orientation of the moments re¬ 
mains almost the same. This phenomenon does not ap¬ 
pear in Fig. [3fa) {N}j = 4). Due to the finite number of 
bath sites used in ED solver, it is not clear whether there 
is indeed a magnetic phase transition in Cl phase. But 
our results indicate that the Cl phase is robust against a 
small change in magnetic order. 

It is interesting to compare the discrepancy between 
the bilayer case and the three-dimensional bulk. To fur¬ 
ther examine magnetic configurations between the bulk 
and the bilayer system, we present the magnetic pat¬ 
tern for ter = —1 in Eig. [5fc). The magnetic moments 
are aligned in close similarity to the AIAO configura¬ 
tion, denoted as F 3 , found in the three-dimensional bulk 
materials The nonzero net moment in the tetrahedron 
is pointing in the [ 111 ] direction due to the failure of can¬ 
celation between moment 1 (in triangular layer) and sum 
of moment 2, 3, and 4 (in the coplanar kagome layer). 
Since the Chern number C = 0, it is still a MI. 


B. Trilayer system 

The ter = — 1 noninteracting band structure for the 
TKT system is shown in Eig. [2][b). Unlike the bilayer 
system, the TKT system has both inversion symmetry 
and TRS. Thus, the noninteracting bands are doubly de¬ 
generate and a metallic ground state is expected with a 
half-filled nearly flat band lying at the Eermi energy. The 
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FIG. 7: Spectral weights along high symmetry lines and local 
density of states for trilayer (TKT) ta = —1. The broad¬ 
ening factor is 77 = 0.02. The black dots show quasiparticle 
spectrum, (a) Cl: U = 0.5; (b) MI: U = 1.0; (c) MI: U = 3.0. 


ground state of the U = 0 system is a trivial conductor 
(denoted “C”). The trilayer phase diagram as a function 
of the interaction U is shown in Fig. [6l 

When the Hubbard interaction reaches U ~ 0.2 (in 
CDMFT), the system breaks TRS and the doubly de¬ 
generate flat bands split into two nearly-flat bands carry¬ 
ing opposite Chern numbers. Therefore, half-filled bands 
give a nontrivial Chern number C = ±1. Meanwhile, 
as can be seen in the spectral weight along high symme¬ 
try line and local density of states in Fig. El^a), one can 
identify a finite gap in this regime. Combining these ob¬ 
servation, it is a nontrivial CL This indicates the single¬ 
particle nature of the interaction-induced Cl appearing in 
TKT. The corresponding magnetic order for TKT is still 
non-colinear and close to the HF results in both bulk^ 
and trilayer— systems, as shown in Fig. [8l^a). 

If the interaction strength is further increased to 
U ^ 0.8, a quadratic band touching appears around 
the F point and the band topology becomes trivial [cf. 
Fig. [Tjb)]. Upon further increasing U, another topologi¬ 
cal quantum phase transition occurs. We confirm this by 
measuring the band gaps around the Fermi level using 
both the topological Hamiltonian and the quasi-particle 
Hamiltonian. By diagonalizing the quasi-particle effec¬ 
tive Hamiltonian defined in Appendix [Bl and calculating 
the corresponding Chern number of the filled bands, one 
finds the same change in the band topology found di¬ 
rectly from the corresponding analysis on the topological 
Hamiltonian determined from the single-particle Green’s 
function. 

For U > 0.8, the TKT system becomes a trivial MI, 
and the magnetic configuration resembles (but is different 
from) the AIAO pattern similar to what is found in the 
bulk. However, for any value of U, the CDMFT results 
do not indicate the presence of the magnetic conductor 
phase [labeled “M” in Fig.[6](b)], as predicted by the HF 
theory. But similarly, the Cl phase in the trilayer system 
is more stable in the CDMFT phase diagram than in the 
HF phase diagram, showing that moderate correlations 
and quantum fluctuations tend to favor to the Cl phase. 



FIG. 8 : Magnetic configuration of TKT within a unit cell, 
(a) ta = —1,U = 0.5; (b) ta = —1,U = 4.0; In (a) and (b), 
the magnetic configuration resembles AIAO but with nonzero 
net magnetization in a tetrahedron. The magnetic moment 1 
and 5 are almost parallel, which enhances the ferromagnetism 
within a unit cell and the net magnetic flux through the [ 111 ] 
plane. 


Finally, we note that while in the bilayer system, the 
Cl phase is “born” out of a parent TI state, the same 
is not true of the Cl in the trilayer: The Cl in the TKT 
system is an example of interactions driving a topological 
state without a parent topological state. Combined, the 
bilayer and trilayer CDMFT results show that the gen¬ 
eral conclusion reached within the the HF theory-that 
a Cl can emerge from both gapped and gapless states 
at the non-interacting level-is robust to the inclusion of 
quantum fluctuations. 


V. SUMMARY 

In summary, we go beyond the single-electron approx¬ 
imation to re-investigate the phase diagram of bilayer 
(KT) and trilayer (TKT) thin film pyrochlore iridates 
in an effective Jeff = 1/2 model with on-site Hubbard 
interaction. By applying the cellular dynamical mean- 
field theory (CDMFT) method with an exact diagonal- 
ization cluster impurity solver, the correlation effects in 
the model Hamiltonian we studied have been fully in¬ 
corporated in the time domain, and partially in the spa¬ 
tial degrees of freedom. Importantly, we show that lo¬ 
cal quantum fluctuations will not destroy the magnetic 
order and interaction-driven topological phases of these 
quasi-two-dimensional systems. The effect from differ¬ 
ent choice of bath levels is summarized in Figl3] and a 
comparison of the CDMFT and Hartree-Fock results are 
shown in Fig l3] and Figlbl An interesting result to emerge 
from this comparison is that in both systems the quan¬ 
tum fluctuations captured in CDMFT (but left out of 
HF theory) tend to stabilize the Chern insulator phase 
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by enlarging the region of Hubbard U values over which 
it occupies the phase diagram. It is worth mentioning 
that we also studied other values of to- and in all cases 
found a similar level of agreement between the CDMFT 
and HF phase diagrams to the phase diagrams shown 
in this manuscript. Thus, we conclude that the HF cal¬ 
culations are reliable for predicting the qualitative fea¬ 
tures of the phase diagrams of systems similar to those 
we study here-namely that there is a non-vanishing spin- 
orbit coupling at the non-interacting level of the Hamil¬ 
tonian. In particular, one does not want to rely on a 
purely interaction-driven spin-orbit coupling to generate 
the Chern insulator phase^^ 

The magnetic order that appears above a critical value 
of U is another important feature of our results. The con¬ 
figurations of local moments under the CDMFT calcula¬ 
tions are non-collinear in both bilayer and trilayer sys¬ 
tems and they highly resemble their counterparts in the 
bulk material with the same tight-binding parameters, 
though they are different in detail because of the lowered 
symmetry of the films. For example, the net magnetic 
moment in the unit cell is non-vanishing due to the quasi- 
two-dimensional nature of the system, which brings ad¬ 
ditional anisotropy in the magnetization. Moreover, our 
CDMFT study indicates that the TI phase exists with an 
interaction strength sufficient to break the time-reversal 
symmetry. This correlated topological phase is adiabat- 


ically connected to the topological band insulator with 
the band dispersion and topology given by the “quasi¬ 
particle effective Hamiltonian” and “topological Hamil¬ 
tonian” , as we define in the text. It is worth pointing out 
that although our numerical results suggest positive out¬ 
comes for topological phases to be hosted in many-body 
interacting systems, the energy scales for these phases to 
be detected in experiments are relatively small-being set 
by the gap value of the quasi-particle bands. The extent 
to which the observable topological transport properties 
survive thermal fluctuations and disorder in interacting 
phases is an important open question for future work. 
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Appendix A: Calculation of Chern number and Z2 
invariant 


In Sec. uni we have briefly introduced the approach 
to calculate the topological invariants using the single¬ 
particle Green’s function. Here we provide more details. 
In a noninteracting topological band insulator, the Chern 
number and Z 2 invariant are evaluated with the following 
TKNN^ and the generic Z 2 formulai^ 


C = 


i-ir = 


/ 


dPk 

27r 


Q:=filled 

\/Det[5(Ai)] 


i=l 


Pf[5(A,)] 


(Al) 

(A2) 


where fij = and a^(k) = 

—i{ua{'k)\-^\ua{^)) are the Berry curvature and 
Berry connection, respectively. Here |i4c,(k)) 

are the noninteracting Hamiltonian eigenstates, 
i.e. iL(k)|i^Q,(k)) = £’Q,(k)|i^Q,(k)). The summa¬ 

tion is over filled band index for Ea(k) < Ep. 
Bmn{K) = (7i^(-k)|T|i4n(k)) is the sewing matrix 

entry on the overlap between Bloch state i4^(—k) and 
the Kramer partner of Bloch state Un(k) at TRIM 
points Ai. In the quasi-two-dimensional thin films, there 
are four TRIM points: F and Mi^ 2 , 3 - An alternative 
approach to calculate the Z 2 index is to see it as an 
obstruction to Stokes’ theorem in the half-BZ, 


P2 = 


[ d^kf 12 — f dk • a mod 2, (A3) 

Jbz JdBZ 


where in the Berry connection for ai(k), a time reversal 
smooth gauge is obtained by choosing the state and the 
corresponding time-reversal partner 


|^,(-k))=T|^,(k)), (A4) 


where T is the time-reversal operator. The topological 
order parameters for interacting topological insulators 
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are derived from topological field theory as 


Appendix B: Quasiparticle Effective Hamiltonian 


Cl 


= 

247r2 J 


cPkcLoe^’^'^ 




dqn dqy 


dqr 


P 2 = 


liup 




1 

120 ^ 


(Fk 
dG-^ 


dq^ 


dq„ dqp 


^dG-^^dG-\ ^ 

G—: —G—— mod n. 


du 


dv 


(A5) 


(A 6 ) 


Note that in the above formulas, one integral is over 
all frequency range. It is shown in Refs. j^l^ that 
the finite-frequency Green’s function G{iuj^ k) is topolog¬ 
ically equivalent to the zero-frequency Green’s function 
G{iuj = 0 ,k). As a consequence, the above formulas can 
be simplified by using the zero-frequency Green’s func¬ 
tion to reduce the integral by one dimension. Further¬ 
more, one can define a “topological Hamiltonian” ^ 

iftop(k) = -G-i(0,k). (A7) 


Since the zero-frequency Green’s function is Hermitian, 
the eigenvalues for the “topological Hamiltonian” are real 


iftop|a,k) =/ia(k)|a,k), (A 8 ) 


In this section, we derive the formalism for the quasi- 



P 2 mod 2 

P 3 mod 2 

ni = 0 

1 

1 

ni = 1 

0 

712 = 0 

1 

1 

712 = 1 

0 

773 = 0 

1 

1 

773 = 1 

0 


TABLE I: The Z 2 index of the 3D bulk pyrochlore iridates at 
U<6.11. 


tion can be written as 

G'a/3(w,k)= (^[(w+//+i(5)-l-t(k)-S(w,k)]“^^c«/3 (Bl) 

where a, (3 are some arbitrary quantum numbers (or¬ 
bital, spin, or sites in a unit cell, etc), and t(k) is the 
Fourier transformed hopping integral. The self-energy is 
in general a complex matrix, but from the Lehman rep¬ 
resentation, we have 

Gt(0,k) = G(0,k), (B2) 


where a is a “band” index. Filled “bands” have eigen¬ 
values /ia(k) < 0. We construct the Berry connection 
and apply Eq. (jA5[) and Eq. (jA3[) with the noninteracting 
Hamiltonian eigenstates |i 4 Q,(k)) replaced by topological 
Hamiltonian eigenstates |(a,k) to evaluate both the first 
Ghern number and Z 2 invariant in pyrochlore oxides thin 
films, which do not have inversion symmetry in general. 
The numerical evaluation of the integral over BZ is based 
on lattice discretization of the gauge field 

We have numerically benchmarked Eq. (jA5[) in the 
three-dimensional bulk pyrochlore iridates by evaluat¬ 
ing the strong Z 2 topological index at U = 0,6,6.11^ 
which has eight time-reversal invariant momentum points 
(TRIM). The reciprocal lattice vector is written as 


and 


St(0,k) = S(0,k). (B3) 

At general complex frequency, we can separate the self¬ 
energy into Hermitian part and anti-Hermitian part as 

S"(^,k) = (S(^,k) + S(V,k)t)/2 (B4) 

S^(^,k) = (S(^,k)-S(V,k)t)/ 2 . (B5) 

The quasiparticle band structure can be defined as^ 


Det 


(cj + /i) • 1 — t(k) — (cj, k) 


= 0 . 


(B 6 ) 


K = nibi-f n2b2 + nsbs, ^ 1 , 2,3 ^ ^ (A9) 

Where bi, b 2 , bs are basis vectors correspond to real lat¬ 
tice vectors ai, a 2 , as. The strong topological index can 
be obtained by P 3 = {P 2 {ni = 0 ) + P 2 (ni = 1 )) mod 2 for 
i = 1, 2, 3, where P 2 is evaluated from Eq. (jA3p . With the 
results summarized in Tabl^ one can verify the strong 
topological insulator index (1;000) in the 3D bulk py¬ 
rochlore iridates. 


To analytically solve Eq. (jB 6 |) around the Fermi level, we 
expand the Hermitian part of the self-energy up to first 
order in cj, 

CJ + ... (B7) 

u ;=0 

Plugging the above into Ea. IBBII . now we have 


S"(w,k) = S"(0,k) + 


aS"(w,k) 


du} 
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Det 


55]^(cj,k) 


duj 


uj =0 


(t(k) -^■l + S"(oj,k)) 


= 0 . 


(B8) 


Ea. dBSIl can be converted into an eigenvalue problem of 
an artificial Hamiltonian. First we diagonalize the matrix 
to the first order of cj in Eg. fjBSp as 


B(k) = 


/ai(k) 0 0 

0 a2(k) 0 

\ 0 0 ■■■ 

as"(w,k) 


(B9) 


f/(k) 1 


duj 


oj=0 


[/t(k). 


For simplicity, we drop the momentum index k in H(k) 
and t/(k) [!/’*'(k)], so that Eg. fjBSp can be rewritten as 


Det 


CJ 


B-U (t(k) - M ■ 1 + S"(0, k)) 


= 0 . 


Let us rewrite 


VF = 


/ ! ^ 
\/^i(k) 

0 


\/^ 2 (k) 


0 \ 
0 


(BIO) 


for the effective quasiparticle Hamiltonian, 

iJ^ff(k) = WU (t(k) + s"(0, k)) U'^W. (B12) 


In the CDMET formalism, since self-energy has no de¬ 
pendence on k, 


E(0,k)«I](0), 


(B13) 


Eg. ()B12p becomes 


^CTMFT(j^) _ -//•! + x;^(o)) u'^w. 

(B14) 

The quasiparticle weight matrix is calculated by applying 
the Cauchy-Riemann equations 



o 

o 



then Eq.(jB8p becomes 



Det 

w • 1 - WU{t{k) - p • 1 + S^(0, k))U'<W 

= 0. 


At this stage, we have defined the “quasiparticle disper¬ 
sion” by solving the eigenvalue problem 

~ ^a/3 


Feff(k)^“(k) =E“(k)V’“(k), 

(Bll) 

where c^o = 


a/3 


(B15) 


duj 


Im 


Cc’= 0 / al3 

S^(u;o)) 


al3 


Re 


(u;o)) 


a/3 


Uo 


LUO 


where cjq = 7r//3 is the first positive Matsubara frequency. 






























